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We study the effect of the interplay between the Rashba and the Dresselhaus spin-orbit couplings 
on the transverse electron focusing in two-dimensional electron gases. Depending on their relative 
magnitude, the presence of both couplings can result in the splitting of the first focusing peak into 
two or three. This splitting has information about the relative value of spin-orbit couplings and 
therefore about the shape of the Fermi surface. More interesting, the presence of the third peak is 
directly related to the tunneling probability ("magnetic breakdown") between orbits corresponding 
to the different sheets of the Fermi surface. In addition, destructive interference effects between paths 
that involve tunneling and those that do not can be observed in the second focusing condition. Such 
electron paths (orbits) could be experimentally detected using current techniques for imaging the 
electron flow opening the possibility to directly observe and characterize the magnetic breakdown 
effect in this system. 
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I. INTRODUCTION 

Transport properties of two-dimensional electron 
(2DEG) and hole gases can be affected in very pecu- 
liar ways by the spin-orbit (SO) coupling. The unusual 
properties of spin transport of these systems are seen as 
promising tools for the development of new spintronic 
devices,^ which would allow us to coherently control and 
manipulate the electrons' spin. This triggered an intense 
activity in the field during the past years. Among the 
SO-related effects in 2DEGs, it is worth mentioning the 
proposal of a spin filtering transistor,i^ the Aharonov- 

Casher oscillation in mesoscopic rings,— and the spin 
Hall effect 6.7.8.9 .10.11.12.13.14.15.16.17.18.19.20.21 

In 2DEGs made either from heterostructures or from 
quantum wells, there are two dominating forms of the SO 
coupling32i The Rashba SO coupling, which arises from 
the asymmetry of the confinement potential of the 2DEG, 
and the Dresselhaus SO coupling, which arises from the 
lack of inversion symmetry of the crystal structure. Both 
types of SO couplings are present in general and their rel- 
ative magnitude depends on the structure and the mate- 
rials used to make the 2DEG. The SO couplings are rela- 
tively weak in AlGaAs-GaAs structures and quite strong 
in In or Sb based semiconductors. There are two impor- 
tant differences between the Rashba and Dresselhaus SO 
couplings. On the one hand, the magnitude of the former 
can be externally controlled by a gate voltage,'^^^ provid- 
ing a interesting new knob to control transport proper- 
ties. On the other hand, the Rashba coupling is isotropic 
while the Dresselhaus coupling depends on the orienta- 
tion of the crystal axes. These two different sources of 
the SO coupling can be experimentally determined using 
different techniques •^i^^'^&i^liS^ 

Both (linear) SO couplings lead to similar electronic 
and transport properties when one dominates. However, 



very interesting effects arise when the two couplings have 
a similar magnitude. In particular, when both couplings 
are equal the spin and the momentum decouple. This 
effect has been proposed as a way to built up a spin tran- 
sistor in disordered systems<^ In addition, it was argued 
that in that case the system contains unusually long-lived 
spin excitations:^^ Indications of the presence of this ex- 
citations have been observed very recently*^ Also, the 
magnetic field anisotropy of the spin relaxation length in 
long wires made from 2DEGs in AlGaAs has been at- 
tributed to the closed values of the two couplingsi^ 



It is therefore interesting to look for new alternatives 
where the effect of the competition between the Rashba 
and Dresselhaus couplings on the transport properties 
can be measured directly. In this work, we analyze the 
effect of such competition on the transverse electron fo- 
cusing signal . ^^'^"^ In Ref. (s^l it was predicted that SO 
couphng leads to the splitting of the odd focusing peaks. 
Since then, this splitting has been observed in differ- 
ent samples2^i2& and discussed by several authors!^2i^i2S 
Here, we show that the splitting of the focusing peaks can 
be used to map out the non-trivial shape of the Fermi sur- 
face of the 2DEGs when both types of linear SO couplings 
are present and have a similar magnitude. In addition, 
we found that the focusing experiment can clearly show 
the tunneling between cyclotron orbits, in direct analogy 
to the magnetic breakdown in bulk materials 



To our knowledge, this is the first example where the 
magnetic breakdown between different cyclotron orbits 
could be directly observed. This could be done by using, 
for instance, the imaging technique developed by West- 
ervelt and coworkers 
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FIG. 1; The Fermi surfaces for /3/a = 0, 0.5, 0.75 and 1 are 
shown in (a), (b), (c) and (d), respectively. The arrows indi- 
cate the spin orientation of the corresponding eigenstate. 



II. SPIN-ORBIT COUPLING IN TWO 
DIMENSIONAL GASES 

A. Bulk eigenstates 

The Hamiltonian of a 2DEG in the presence of both 
Rashba and Dresselhaus SO couphngs is given by 

a (3 

H = 2^ + '^(Pv^^ -PxCTy) + -^{Px'CTx' -Py'CTy') (1) 

where p = {px,Py) is the momentum operator, a and 
/3 are the Rashba and Dresselhaus coupling parameters, 
respectively, and {ci} are the Pauli matrices. The axes 
x' and y' correspond to crystallographic directions while 
X and y are arbitrary directions chosen in a convenient 
way — note that the Rashba term is isotropic and there- 
fore independent of the axes choice. 

Hamiltonian ([1]) can be easily diagonalized by propos- 
ing a solution in the form of a plane wave. The eigen- 
functions and eigenvalues are 

*±(r) = -Le"^'-f^]) (2) 

e(fc) = ^±^{akx+f3ky)^ + {aky+f3kx)^ . (3) 

Here, tan^ = —{akx + (3ky)/{aky + (ikx), A is the sys- 
tem's area and we choose x = x' and y = y' . Figure [T] 
shows the corresponding Fermi surface for different values 
of the ratio (3 /a. The arrows indicate the spin orienta- 
tions of the eigenstates. 

The competition between the Rashba and Dresselhaus 
SOs originates a deviation of the Fermi surface from the 
circular shape. For a = ±/3 the Fermi surfaces recover 
a circular shape shifted from the F point (there is a per- 
fect nesting between the two surfaces) and the spin ori- 
entation becomes independent of k — nevertheless, this 
case has interesting spin properties i^i^Sii^i^ For \f3/a\ ^ 1 
the two different Fermi surfaces do not cross each other. 



The minimum and maximum distances between them in 
fc-space are given by AA:± = ^p-la ± /3|. As we will 
show, these properties have important consequences on 
the transverse focusing signal . 

When an external magnetic field is applied perpen- 
dicular to the sample. Landau levels are formed. In 
that case, a Zeeman term must be included in Hamil- 
tonian ([T]) and p should be replaced by p-f (e/c)A, with 
A as the vector potential. A closed analytical solution 
for arbitrary values of a and /3 is not known (see how- 
ever Refs. pglfeo] ) . For /3 = 0, however, a straight- 
forward calculation^^ shows that the energy spectrum is 
given by = hujcn =F (^^o + (Q^/^c)^2ri) 2 , where n > I, 
LUc — e\ B\ /■m*c is the cyclotron frequency, Ic — (h/muJc)'^ 
is the magnetic length, and EQ = hWc/2 — gjisBz/'^. is the 
energy of the ground state (n = 0). In the limit of strong 
Rashba coupling or large n, (a/^c)\/2n 3> -Eq, the spin 
of the eigenstates lies in the plane of the 2DEG. These 
eigenstates have a cyclotron radius given by 

rl ~ 2nll (4) 

Then, states with different n, and consequently different 
cyclotron radii, coexist within the same energy window.— 
In fact, it is easy to verify that the difference between the 
two cyclotron orbits is 

2a hAk± „ , , 
hujc mrujc 

Equivalent results are found in a semiclassical treatment 
of the problemj^i21i2^i^ 

Eventually, if both a and (3 were non-zero, one might 
hope to be able to gather information on the shape of the 
Fermi surface by measuring Aj'c. As we show below, this 
is indeed the case when transverse electron focusing is 
used to map out Arc as a function of crystal orientation 
or the strength of the Rashba coupling. 

B. Transverse electron focusing 

The usual geometry for transverse focusing experi- 
ments consists of two quantum point contacts (QPCs) 
at a distance L, which are coupled to the same edge of 
a 2DEG (see Fig. H]). Electrons emitted from QPC / 
(injector) are focalized onto the QPC D (detector) by 
the action of an external magnetic field perpendicular to 
the 2DEG. In a classical picture, the electrons ejected 
from the injector are forced to follow circular orbits due 
to the Lorentz force. If the applied magnetic field has 
some arbitrary value, the electrons miss the detector and 
simply follow skipping orbits against the edge of the sam- 
ple. However, for some particular values of the external 
field (Sn), such that the distance between the QPCs is 
n times the diameter of the cyclotron orbit, with n as an 
integer number, the electrons reach the detector. In such 
a case, there is a charge accumulation in the detector 
that generates a voltage difference across QPC D. This 
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FIG. 2: Scheme of the focusing experiment in the presence 
of SO coupling (we take (3 = for simplicity). An electron 
injected at QPC / can follow one of two different orbits (Oi 
and O2) depending on its spin orientations. These two or- 
bits lead to the splitting of the first focusing peak since the 
corresponding focusing fields and _Bi,2 are different. Af- 
ter bouncing off the edge the electron that followed orbit the 
smaller (bigger) orbit continues in the bigger one (smaller one) 
so in both cases the electron arrives to the detector D for the 
same focusing field B2, leading to a single peak. 



gives voltage peaks as the external field is swept through 
the focusing fields B„^iSi^ 

In a quantum-mechanical description, the scattering 
states in the two QPC are coupled by the Landau levels 
of the 2DEG. As the Landau eigenstates have a char- 
acteristic length (the cyclotron radius Tc) that depends 
on the applied field, there is also a matching condi- 
tion for B = Bn- In this context, the main features 
of the magnetic-field dependence of the measured sig- 
nal are contained in the transmission T between the two 
QP(-i33 — typical experimental setups include also one or 
two Ohmic contacts at the bulk of the 2DEG which are 
used to inject currents and measure voltages<^ 

As shown in Refs. [131 and [35], in systems with ei- 
ther Rashba or Dresselhaus spin-orbit coupling (but not 
both) the first focusing peak splits in two. Such splitting, 
for /3 = 0, is given by Eq. ([5]) or, in terms of the mag- 
netic field, by Ai? — 'iam*c/heL, which is independent 
of B. Furthermore, each peak corresponds to a differ- 
ent spin projection of the electron leaving the emitter 1^ 
Once again this can be understood using a classical pic- 
ture plus the fact that there are two Fermi surfaces, even 
though the semiclassics is not triviali^^i^i^i^iSi^i^ 

This simple mechanism, which is able to spatially sep- 
arate the two spin orientations of a electron beam, was 
recently used^ to study the current's spin polarization 
associated with the "0.7" anomaly in QPCs and it was 
also suggested^ as a tool to study spin polarization of 
the flowing current in adiabatic QPCs due to SOi^ 



C. Numerical solution 

As mentioned above, we are interested in calculat- 
ing the conductance between the two lateral QPCs. In 
the zero temperature limit this conductance is just e^/ft- 



times the transmission coefficient T between the two con- 
tacts evaluated at the Fermi energy. Even in the absence 
of the spin-orbit interaction, it is not possible to obtain 
an analytical solution of the problem when there is an ap- 
plied perpendicular magnetic field. Hence, we calculate 
T[Ep) numerically using a discretized system ("tight- 
binding-like" model) where the leads or contacts can be 
easily attached. The Hamiltonian of the system can then 
be written as H = Hq + Hsoi where 



(n,m) ,(T 



-h.c. 



(6) 



Here, cj^^ creates an electron at site n with spin a (t or 
I in the z direction) and energy = At — agjiBBz/'i, 
t — h? /2m*a^, and ao is the effective lattice param- 
eter which is always chosen to be small compared to 
the Fermi wavelength. The summation is made on 
a square lattice, where the position of the site n is 
UxX + riyy, where x and y are unit vectors in the x 
and y directions, respectively. The hopping matrix ele- 
ment t„„ is nonzero only for nearest-neighbor sites and 
includes the effect of the diamagnetic coupling through 
the Peierls substitution's. For the choice of the Landau 
gauge =iexp (-in^27r0/(?!)o) and *„(„+§) =t with 

(j> = GqB the magnetic flux per plaquete and 0o = hc/e 
the flux quantum. 

The second term of the Hamiltonian describes the spin- 
orbit coupling, 

n 

where = a/2ao, Xd = /3/2ao, A^; = {XR + iXne^'^'^) 
and Ay = — (Aj^e^'^'^' + iAfl) and (p is the angle between 
the crystallographic axis x' and the x axis (normal to 
the edge of the 2DEG). In the second term the Peierls 
substitution is made explicit. 

Each lateral contact is described by a narrow stripe 
with a width of Nq sites and, for simplicity, no spin-orbit 
coupling. They represent point contacts gated to have 
a single active channel with a conductance 2e^/h (for 
details see Ref. 0|). To obtain the conductance between 
the two contacts we calculate the retarded (advanced) 
Green's function matrix, Q^'-°-K Because of the lift of the 
spin degeneracy the Green's function between two sites i 
and J has four components Qia,ja' ■ 

The zero-temperature conductance is then 
obtained using the Landauer formula, 6*12 — 
(eV/i)Tr[r(i)^;'T(2)^'^], evaluated at the Fermi en- 
ergy. Here F^^'' =i[S^ — S^] is the "coupling matrix" to 

the contact N and S^"'' the corresponding self-energies 
of the retarded (advanced) propagator. 
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FIG. 3: (a) Splitting of the first focusing peak as a function 
the crystallographic angle ip for Ef = 23meV, m* = 0.055mo, 
L = 1.5/im, a = 15meVnm, and /3/a = 1/3. The dashed lines 
correspond to Eq. @ and B/ = 165mT; (b) and (c) show the 
focusing peaks for the minimum and maximum splitting; (d) 
scheme of the Fermi surface indicating the difference Ak{if). 



III. SPLITTING OF THE FOCUSING PEAKS 



and a scale factor l'^. Since in our case there are two 
different Fermi surfaces, there are two real-space orbits 
whose radius difference is then given by Ar{ip) = llAk{ip) 
where Afc(i^) = {2ni* / fi^) y/ a'^ + f3'^ - 2af3sm2ip is the 
orientation-dependent difference of the two wave vectors 
of the two Fermi surfaces (see Fig [3]). Here, in deter- 
mining the cyclotron radius, we have neglected the fact 
that the velocity is not parallel to k (normal injection 
does not always correspond to ky — 0). In the limit 
Ak{(p)/kp <C 1, as it is the case here, this is an excellent 
approximation. 

Then, the peak position is given by B±{ip) = B f{Lp) ± 
AB{ip)/2, where Bf{(p) 



^V2Em* and 



ABiip) 




- 2 



sin 2^3 (9) 



is the magnetic-field splitting of the first focusing peak. 
This is indicated in Fig. [3] with dashed lines. 

The agreement with the exact numerical result is ex- 
cellent, showing that the simplified semiclassical picture 
is very accurate in this regime. From Eq. ([9]), we see 
that a measure of the peak splitting is a direct way to 
measure a and /3. One possible way to do it would be to 
use different sets of pairs of QPCs oriented in different 
angles with respect to the crystallographic axis. 



A. Dependence with the crystal orientation 



B. Additional peak: "magnetic breakdown" 



Figure [3] shows the splitting of the first focusing peak 
as a function of the crystal orientation (defined by the 
angle (p) with respect to the edge of the sample. The 
splitting shows a simple oscillatory behavior, whose an- 
gle dependence can be fully understood in terms of the 
shape of the Fermi surface and a simple semiclassical ar- 
gument presented below. It is worth mentioning that the 
semiclassical description of the orbits is far from trivial in 
the presence of spin-orbit coupling ..^L^i^i^i^i^i^iS In 
particular, in the presence of both Rashba and Dressel- 
haus couplings, there is an extra complication due to the 
possibility to have mode conversion points (points where 
the spin-orbit field cancels) .^i^i^iSi^i^ This will be 
important for the effect discussed in the following sec- 
tion. 

Here, however, we can argue that the SO coupling is 
sufficiently strong so that the spin follows the momentum 
adiabatically. We can then use the usual semiclassical 
description for a band structure given by Eq. ([3]). In 
that case, the semiclassical equations of motion are given 
by 



r = t> = iVfce(fc) , 
n 



k — —v X B . 

he 



(8) 



For B = Bz, the solution of these equations implies 
that k{t) moves along the Fermi surface while r{t) — 
r{0) + z X {k{t)~k{0)) ll and, as usual, the real-space 
orbit is related to the one in fc-space by a 7r/2 rotation 



Let us now consider a different situation where the 
crystallographic orientation is kept fixed (we take x' = 
X so that Lp — Q) but the magnitude of the spin-orbit 
coupling (a or (3) is changed. This is shown in Fig. |4l 
When one of the SO couplings dominates the two peaks 
structure is clearly seen. Consider the case of small a/ (3 
in Fig. m as a/ (3 increases the splitting also increases in 
agreement with Eq. However, for a/ (3 ~ 1, a third 
peak develops at B = Bj. The amplitude of this peak 
increases at the expenses of the other two and becomes 
the only peak for a = [3. The fact that there is only 
one peak when a — (3 is quite clear from the fact that in 
that case the two Fermi surfaces are circular and have 
the same radius. 

The transition from two to three peaks can be under- 
stood in terms of tunneling between cyclotron orbits, in 
the same spirit as the magnetic breakdown between band 
orbits in bulk metals3 For a ~ /3, the "gap" in fc-space 
between the two Fermi surfaces, that is the minimum dis- 
tance between them Afc_ ~2m*\a — (3\/'h? , is very small 
and the magnetic field can induce a tunneling transition 
between both orbits. Therefore, if f{B — B*) describes 
a single focusing peak centered around B* , the complete 
focusing signal is expected to behave as 

{.f{B - B+) + f(B - B^)) il-p) + 2pf{B-Bf) , (10) 

where p is the tunneling probability, which can be esti- 
mated using a Landau-Zener-type argument. Following 
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FIG. 4: Magnetic breakdown between different cyclotron or- 
bits, (a) Evolution of the first focusing peak as a function 
of a and /3 when fixing the greater in 15 meVnm. The split- 
ting evolves as expected from Eq. ((9)1 until the two coupling 
have similar values. Close to this point a third peak appears 
as a consequence of the quantum tunneling between the two 
orbits. For q = /9 we have only one orbit, AA: = 0, and then 
only one peak; (b), (c) and (d) show the focusing peaks for 
different values of the SO coupling parameters; (e) focusing 
signal as a function of (p for a//3 — 0.75; (f) fitting of the 
amplitude of the central peak with Eq. (fT2)l . 



Refs. [mil, a rough estimate for p is 



p = exp 



-nli 



Afc3 



(11) 



where 1/a and 1/6 are the curvature radii of the two 
Fermi surfaces in the tunneling region, ie. close to the 
minimum gap point. For a/P ^ 1, this can be approxi- 
mated by 



p — exp 



Tnn*kFll (a- Pf 



exp 



-7- 



(12) 

with X — a/ (3 and 7 = ■nm*kpl'^[3 /h^ . This expression 
fits the numerical data [see Fig. [3f] up-to a factor of 0.6 
in 7. Notice that the fitting function is not a Gaussian 
and that p is independent of the crystallographic angle 




FIG. 5: Imaging of the cyclotron orbits (see text), (a) and 
(b) show the total conductance Ginj,tip for two different crys- 
tallographic orientation, tp = 0.757r and 0.7257r, respectively, 
13/a = 0.75 and L = 1.5/im. Note that in both cases there 
are three peaks at the first bounce against the edge of the 
sample while in the second bounce the central peak is missing 
(due to interference) only in (a); (c) and (d) show the spin 
resolved conductances, Ginji.tipT ~ Ginji.tipi- The absence of 
a sign change along the orbit is a signature of the magnetic 
breakdown. 



(/?; it only depends on the local properties of the Fermi 
surface around the minimum gap. 

At this point it is worth mentioning that the mag- 
netic breakdown of the cyclotron orbits has previously 
been used to explain the anomalous behavior of the 
magnetoresistance oscillations in systems with spin-orbit 
coupling.'*^ This interpretation has been challenged very 
recently^^i^i^ arguing that the dynamics of the spin can- 
not follow, in that case, the momentum. Our case, how- 
ever, is different as we are in the strong spin-orbit limit 
and the magnetic breakdown interpretation is appropri- 
ate. 

In order to explicitly show the magnetic breakdown 
of the cyclotron orbits, we plot in Figs. [5^ and [SJs the 
conductance Ginj,tip from the injector to a conducting tip, 
located above the 2DEG, as a function of the tip position 
(see Ref. 3J| for details) and for two different crystallo- 
graphic orientations, Lp = O.TStt and 0.7257r, respectively. 
This essentially corresponds to calculate the probability 
for an electron injected through QPC I to reach a given 
point in the 2DEG and hence brings information about 
the orbits followed by the electrons. The images obtained 
in this way are similar, although with a much better reso- 
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FIG. 6: Focusing signal at the second focusing peak: (a) as 
a function of a (/3) for /3(a) = ISmeVnm and i^j = 0. As 
expected there is only a single peak, except for a//3 ~ 0.75 
where there is a destructive interference between the differ- 
ent orbits shown in Fig. [T] Notice this point correspond to 
the case where there are three peaks at the first focusing con- 
dition; (b) as a function of the crystallographic angle ip for 
a//3 ~ 0.75. 



lution, that the ones we would have obtained by simulat- 
ing the presence of a tip as a scatterer (th e experimental 



technique developed in Refs. 

Although it is difficult to distinguish the orbits along 
the full path, three different orbits are apparent in both 
cases close to the first bouncing point, which correspond 
to the first focusing condition. As each of these or- 
bits have a spin projection associated with it, we plot 
in Figs. [S}: and [SJi the difference between the spin re- 
solved conductances, Ginjt.tipT^Ginjt.tipi, wheic the spin- 
quantization axis corresponds to y. This allows us to 
follow the direct orbits and the ones that involve tun- 
neling. For this, it is important to take into account 
that for /3 = 0, we would observe a change of sign of 
Ginjt,tipT ~ Ginjt,tipi in the middle of the orbit as the 
spin rotates from | t) to I i)- Here, the fact that there 
is an orbit in which this quantity does not change sign 
is an indication of the tunneling from one orbit to the 
other. To support this interpretation, we also show the 
semiclassical orbits (dashed lines) expected for an elec- 
tron injected with a given spin polarization along the y 
axis. To include the orbits that involve tunneling, we 
simply change from one orbit to the other at the position 
related to the minimum gap in fc-space (see Fig. [7]). 

An interesting effect occurs at the second focusing peak 
(second bounce in Fig. O. In such a case, the peak struc- 
ture results from the interference of several paths and 
hence destructive interference can results for particular 
values of the parameters. In particular, we see in Fig. [5] 
that for ip — O.TStt the central peak is missing. In order 
to understanding the origin of this effect, we show in Fig. 

the focusing signal as a function of a//3 for 1^9 = and 
the same microscopic parameters as in Fig. For almost 
all values of a//3 there is a single peak (as in the /3 = 
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FIG. 7: Semiclassical orbits (in real space) containing two 
bounces for an electron injected at the point indicated by the 
arrow with its spin polarized along the edge of sample, (a) 
Direct orbit (no tunneling) and (b) orbit with two tunneling 
events, indicated by the circles. Here <^ defines the angle 
where tunneling occurs; (c) superposition of orbits (a) and 
(b). Notice that after the second tunneling event the two 
orbits overlap; (d) and (e) orbits that involve one tunneling 
event; (f) superposition of all previous orbits. The letters 
A,B,C and D identify the different locations of the bounces. 



case). However, for a//3 ~ 0.75, when the first focusing 
peak shows a three peak structure (see Fig. |4]), the cen- 
tral peak disappears while two satellites peaks emerge. 
When setting a//3 = 0.75 this effect has an oscillatory 
behavior as a function of the crystal orientation. This is 
shown in Fig. [6)3 (see also Fig. Ilje) for a comparison 
with the first focusing peak structure). 

As mentioned above, the origin of this modulation of 
the amplitude of the second focusing peak is the inter- 
ference between the different orbits that contribute to 
the signal. These orbits, for an electron injected with 
its spin pointing along the y-axis, are shown in Fig. [T] 
The regions in real space where tunneling between the 
two orbits occurs are indicated with circles. The posi- 
tion of these regions depends on the crystal orientation. 
Figures [7^ and [7)3 show the orbits that contribute to the 
central peak of the second focusing peak, while Fig. [71; 
shows the overlap of the two orbits. Figures [Tli and [7^ 
show the orbits that contribute to one of the satellites 
in the second focusing peak. The letters A,B,C and D 
identify the different locations of the bounces and then 
the different focusing peaks. 

As the phase acquired due to the tunneling between 
orbits is independent of in order to account for the 
angle dependence of the interference pattern we need to 
include the orbital phase acquired by the electron along 
the different paths. In addition, we notice that after the 
second tunneling event the two orbits in Fig. [Tj: follow 
exactly the same path. Since they accumulate the same 
orbital phase from there to the next bouncing point, if 
they interfere destructively right after the tunneling, they 
will do it along the final part of the orbit. This is clearly 
seen in Fig. [5]where for = 0.757r the orbit correspond- 
ing to the central peak in the second focusing peak has 
disappeared completely after the second tunneling event. 
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The central peak amplitude is then given by 

A = I ^/P'd + /^e'(2e'+»-^) |2 (13) 

where Pd and Pt are the probability of the direct and 
tunneling paths, respectively [paths (a) and (b) in Fig. 
[7], dt is the phase acquired by the electron due to each 
tunneling event (we assume them to be equal), and 9 orb — 
6i — 02 is the difference of the orbital phases of the two 
paths. In the semiclassical picture, the orbital phase of 
each path is given by 1/hJ^p-dr, where the integral 
is done along the classical path C, p is the canonical 
momentum and r the position vector satisfying 



dH 
dp 



P 



dH 
dr 



(14) 



In the strong spin-orbit limit, the Hamiltonian H is given 

1j ,,38.57.59 



p2 



H = ^± ^{aP,'+PPy'? + {aPy,+l3P,.y (15) 

with P = p + {e/c)A. Using a symmetric gauge, 
it is straightforward to show that 1/h^^p ■ dr — 
(1/2Z^) J rpd(/), where rc is the radius of the correspond- 
ing cyclotron orbit. The orbital phase difference between 
the two paths is then 



rb = 





(16) 



with if) as the angle where tunneling occurs (see Fig. [7]) 
and Ti, i = 1,2, as the radii of the orbits. This integral 
can be calculated analytically for arbitrary values of a 
and (3 in terms of elliptic integrals. However, for a/ (3 :^ 1, 
it can be approximated by 



ft2 



■cos((p- -). 



(17) 



Figure [8] shows the total central peak amplitude ob- 
tained numerically as a function of ip. Since for Lp it /A 
the satellite peaks merge into the central peak, in Fig. 
[8)3 we plot the central peak amplitude subtracting the 
satellite contributions. This can be done because the 
orbits that contribute to the central peak for all angles 
are those where the electron arrives with the same spin 
orientation it had at the injector. A fitting of the peak 
amplitude, using Eq. (fT5|) with Pd and Pt as the only 
fitting parameters, is shown in Fig. [8h. The dotted line 
was obtained using the approximated expression for Oorb 
[Eq. (|17p ]. while the dashed line corresponds to the exact 
expression [Eq. (|16p ]. Once again, the agreement is very 
good, given further support to the magnetic breakdown 
interpretation. 

It is easy to verify that for the case of the satellite 
peaks, Oorb takes the same value than in the case of the 
central peak. However, as in this case both interfering 
paths contain a tunneling event, 9t drops out. Then, 




FIG. 8: (a) Amplitude of the second focusing peak (indicated 
by a dashed line in Fig. [Gjl as a function of the crystallographic 
orientation. The oscillation is due to quantum interference 
between different paths (see Fig. [Tjl; (b) same as above but 
without taking into account the contribution from the satellite 
peaks. The dashed (dotted) line corresponds to a fitting using 
the exact (approximated) value of Oorb, given by Eq. (|16|l [Eg. 



m] 



the fact that there is destructive interference in the cen- 
tral peak when the satellites have there maximum values 
indicates that 9t — 7r/2 as expected from the magnetic 
breakdown picture. In the above analysis, we did not 
consider explicitly the phase acquired by the spin degree 
of freedom. In general, such phase has a nontrivial depen- 
dence with the geometry of the orbit.— Here, 
since we assumed that the spin follows the SO field, that 
is, it rotates around the z axis, this phase is not relevant. 
However, special care should be taken when considering 
the effect of the tunneling between orbits as the spin may 
have an additional rotation. As this case involves a mode 
conversion point, a semiclassical analysis of this phase is 
not simple — we are not aware of a simple way to esti- 
mate it. Therefore, the fact that our results indicate that 
the phase due to tunneling is just the usual 7r/2, is an 
indication that the spin phase is or tt. Further work is 
needed to clarify this point. 



IV. SUMMARY 

We have shown that the interplay between the Rashba 
and Dresselhaus couplings introduces new effects on the 
transverse electron focusing. The most interesting aspect 
is the appearance of additional structure of the focusing 
peaks related to the magnetic breakdown of the cyclotron 
orbits when the two SO couplings have similar magni- 
tude: the two sheets of the Fermi surface lead to different 
paths in real space and, as we have shown, the tunneling 
between different paths generates new structure in the 
focusing peaks. In addition, interference effects between 
these paths lead to an oscillatory behavior of the second 
focusing peak amplitude as a function of the orientation 
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of the crystallographic axes. We have shown that the ob- 
served interference effects are dominated by the orbital 
phase accumulated along the different paths. This is so 
because in this regime, the spin adiabatically follows the 
momentum and the associated Berry phase cancels out 
while the spin phase due to tunneling seems to be irrele- 
vant here. One could however envision a different regime 
where the spin dynamics would be important for which 
a deeper understanding of the spin dynamics during the 
tunneling as well as of the semiclassical description of the 
problem is needed. 

The magnetic breakdown of the orbits as well as the 
interference effect could be directly observed using some 
of the techniques recently developed for imaging the elec- 
tron flow.— i^i^i^ In this work we have been mostly con- 
cerned with the dependence of the different effects on the 
crystal orientation. In practice, the experimental obser- 
vation is a bit cumbersome as it requires us to tailor 
different QPC setups in different orientations. An alter- 
native to this could be the use of in-plane magnetic fields 
to modulate AA:_ as well as and k^. For instance, for 



a ~ /3, an in-plane field applied along the x' — y' direc- 
tion will essentially control the magnitude of the gap in 
fc-space, and then the probability for tunneling. On the 
other hand, in-plane field along the x' + y' direction can 
be used to modulate /c+ and fc_ and then the interference 
pattern. This also have the advantage of keeping fc+ + k- 
constant and then the focusing condition. 

Finally, it is interesting to note that while the splitting 
of the focusing peaks allows for the spatial separation of 
the spin components as in a Stern-Gerlach device,— 
the presence of the magnetic breakdown between orbits 
provides a way to separate a given spin component in a 
superposition of two spatially separated orbits. 
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